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Abstract 

With an effective Lagrangian approach, we analyze several A^A^ — ?> NNtttt channels by including 
various resonances with mass up to 1.72 GeV. For the channels with the pion pair of isospin zero, we 
confirm the dominance of A^*(1440) — )• Na in the near threshold region. At higher energies and for 
channels with the final pion pair of isospin one, we find large contributions from A^*(1440) — )• Avr, 
double-A, A(1600) Af*(1440)7r, A(1600) An and A(1620) Avr. There are also sizeable 
contributions from A — )• Avr, A — ?> A^vr, A^ — ?> Att and nucleon pole at energies close to the 
threshold. We well reproduce the total cross sections up to beam energies of 2.2 GeV except for 
the pp — 7> ppn^TT^ channel at energies around 1.1 GeV and our results agree with the existing data 
of differential cross sections of pp pp-K^ir^ , pp nnir^-K^ and pp — ?• ppir^iT^ which are measured 
at CELSIUS and COSY. 
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I. INTRODUCTION 



Double pion production in both pion- and photo-induced reactions has been an intriguing 
field to study baryon spectrum and given insight to the properties of strong interaction jl!,!^. 
These reactions close to threshold are also an interesting area to test chiral symmetry and 
have been extensively explored experimentally [sl and theoretically j^. Recently the double 
pion production in the electro-production off protons has advanced a important step [5]. All 
essential contributions are identified from the data and the major isobar channels are well 
determined. On the other hand, as another fascinating platform for studying resonances 
properties, double pion production in nucleon-nucleon collisions has been accurately mea- 
sured at the facilities of CELSIUS and COSY in the past few years, and the comprehensive 
data of various differential cross sections are obtained up to beam energies 1.3 GeV jsHl. 
However, on the theoretical side, the study of this reaction is scarce. The state-of-art one is 



still the Valencia model calculation 



12| of more than ten years ago after some much earlier 



calculations of one pion exchange (OPE) model 13|] of more than 45 years ago. Thus a more 
comprehensive analysis matching the modern data is very necessary. 

The early OPE model, which mainly focused on the old data at beam energies of 2.0 GeV 



and 2.85 GeV 



14j . included two types of diagrams with the final two pions produced from 



a single and two baryon line(s), respectively. It used the amplitudes of tcN — t- ttttN and 
irN — )■ ttN extracted from limited data and the off-shell corrections were considered under 
several assumptions. It did not account for the explicit production mechanisms of double 
pion and other exchanged meson besides vr-meson. The Valencia model is characteristic by 
the the dominance of A^*(1440) — > Na in the near threshold region in the isospin allowed 
channels while the double- A and A^*(1440) — > Avr rise up at higher energies and in channels 
where iV*(1440) — )■ Na is forbidden by the isospin conservation. Recently, the experimental 



data 



i, |9|] confirm the predicted behavior close to threshold and this makes pp — > ppvT+vr" 
and pp ppTc^TT^ good places to study A^*(1440) whose structure is still controversial. 
Current data seem to show a weaker A^*(1440) — > Avr than that listed in Particle Data 
Group [l5| and support the explanation of A^*(1440) as the monopole excitation of the 
nucleon. Contrarily, the case is much more complicated at higher energies. The dominance 
of the double-A mechanism in the Valencia model results in that the total cross section of 
pp — )■ ppiT^TT^ is about a factor of four larger than that of pp — > nrm^-K^ , while the new 
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exclusive and the old bubble-chamber data are consistent to conclude an approximate equal 
value of these two channels. The isospin decomposition unambiguously reveal that more 
isospin 3/2 resonances besides A is required to explain the data p, |8|], and this is also the 
reason that the Valencia model including simply the A^*(1440) and A achieved merely a 
rough agreement in most channels. Indeed, at higher energies the contribution from higher 
lying resonances, especially those having large double pion decay channels, should become 
relevant. The recent detailed measurements performed by CELSIUS and COSY make the 
further exploration of these problems possible. 

In the present work, we try to incorporate the resonances with mass up to 1.72 GeV in an 
effective Lagrangian model with the motivation to give a reasonable explanation to the six 
isospin channels of A^A^ — )■ NNtttt simultaneously and get better understanding of dynamics 
for this kind of reactions. Our paper is organized as follows. In Sect. [Tll we present the 
formalism and ingredients in our computation. The numerical results and discussion are 
demonstrated in Sect. IIIII and a brief summary is given in Sect. 



II. FORMALISM AND INGREDIENTS 

We consider nearly all possible Feynman diagrams as depicted in Fig. 1, and exchanged 
diagrams are also included. We use the commonly used interaction Lagrangians for irNN, 
ttAA, r]NN, aNN and pNN couplings, 

^ fnNN- 



-nNN 



-iV757Mr-9^vfiV, (1) 



C.AA = ^ A'^757^r • 9^7fA, + h.c, (2) 

^r,NN = -igr,NNN-f5r]N, (3) 

-CaNN = gaNwNaN, (4) 

CpNN = -gpNNNin^ + 7r^fT^.9")f • ^N. (5) 

At each vertex a relevant off-shell form factor is used. In our caculation, we take the same 
fonn facto, a. that used the Bon. potential .odelH. 



Fm [kM) = [ .2 _p ) ' (6) 
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FIG. 1: Feynman diagrams for A'^A'" NNtttt. The solid, dashed and dotted hnes stand for the 
nucleon, mesons and intermediate a{or p)-meson. The shading histograms represent the interme- 
diate resonances or nucleon poles. In the text, we use R — ?■ NM, Rl — >■ R2M and double-i? to 
label (1)(2), (3) (4) and (5) (6), respectively. 



with n=l for vr- and r^-meson and n=2 for p- meson, fcjvf, and are the 4-momentum, 
mass and cut-off parameters for the exchanged meson, respectively. The coupling con- 



stants and the cutoff parameters are taken as jl6l-ll8l]: 11^^/^'^ — 0.078, g'^j^jy/in = 0.4, 



sInn/'^^ = 5.69, glj,j^/^7c = 0.9, A. = A, = 1.0 GeV, A, = 1.3 GeV, A, = 1.6 GeV, and 
K = 6.1. We use Jt^aa = ^f-nNNl'^ from the quark model [l|, isl. The mass and width of 
(T-meson are adopted as 550 MeV and 500 MeV, respectively. 

We include all A^* and A* resonances with spin-parity 1/2^, 3/2^, 5/2^ and mass up 
to 1.72 GeV listed in Particle Data Group (PDG) tables 15|. The resonances with further 
higher masses are expected to give negligible contributions in the energy region considered 



here and their two pion branching ratios have large uncertainties, so we do not inc 
at present. The effective Lagrangians for the relevant resonance couplings are 



19 



ude them 



2Q|, 



^llNR = QrjNRN-f^ljR + h.c, 
J-'aNR = QaNRNaR + /i.C, 



(7) 

(8) 
(9) 



'^pNR = 9pNrN^^t ■ p^^R + h.c, (10) 

^HIr = g.ARA^f ■ d^m + h.c, (11) 

JO^J^r^ 9nNRNr -TrR + h.c, (12) 

^vNR = gnNRNrjR + h.c, (13) 

J^Inr = 9pNrNi,{i, - ^)r • pi^R + /i.e., (14) 

^i-Ail = dnAR'^p.-y^r ■ 9^7fi? + /i.e., (15) 

= ^.iVfliVf • d^TTRf, + h.c, (16) 

^Inr = gvNRNd'^vR^. + h.c, (17) 

^piVR = gpNRN'y^f ■ p^'R^ + /i.e., (18) 

^^^A^ = gnAR^'^^bT ■ nRi^ + h.c, (19) 

45^(1440)ii = ^^iV*K^^T • a^Tfi?^ + /i.e. , (20) 

= g.NRN-f5j^r ■ d^d^nR, + h.c, (21) 

^Inr = gpNRNf ■ pi^R^ + h.c, (22) 

^IaR = g^KAR^^lvT ■ d^T^R^ + h.c, (23) 

^Inr = 5^iViiA^757M^ • d^d-'d^nR^x + h.c, (24) 

= ^P^^AT (^p^;, - ^^^^^) r ■ p^^i?,. + /I.e., (25) 

^ITr = gaNRNd'^d'^aR^, + h.c, (26) 

i'i = ^.tAr • d'^d-'d'^nR^x + h.c, (27) 

= ^.iViiA^T ■ ^^^'^Tfi?^. + h.c, (28) 

>C'aR = ^.rARA^75T • d^d^d'^TfR^x + h.c, (29) 

For the Resonance-Nucleon-Meson vertices, form factors with the following form are used: 

^■r(*L)=(^^)". (30) 
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with n=l for A^* resonances and n=2 for A resonances. We employ A* = A* = A* = A* 



= 1.0 for all resonances except A* = 0.8 for A* 
factors B{QN*An) in the A^*(1440)-A-7r vertices 



600). We also use Blatt-Weisskopf barrier 



2l|, 



B{Q 



N*Att 



(31) 



Qn'-Att + Qo 

Here Qo is the hadron scale parameter Qo = 0.197327/-R GeV/c, where R is the radius of 
the centrifugal barrier in the unit of fm and is tuned to be 1.5 fm to fit the data. Qn* An 
and Pn* An is defined as, 

n2 _ (-^A^ + Sa- -Stt) fr,r,\ 

'^N*An — ~ ■5a, {OA) 



Ac* 



m: 



2\2 



N*An 



4:m%, 



— m 



A? 



(33) 



with Sx being the invariant energy squared of x particle. Because the mass of a-meson is 
near the two-vr threshold, the following Lagrangians and form factor are employed for the 



cr-TT-TT vertex 



pTTTT 



9pnnT^ X a^TT ■ p^. 



(34) 
(35) 

(36) 



A2+g2 

where q is the relative momentum of the emitted vr-mesons. We use A = 0.8 GeV and 
Aq = 0.12 GeV^ to normalize this form factor to unity when vr- and a- meson are all on- 



shell. The decay width of a — )■ vrvr and p — )■ vrvr yield g^^^ = 6.06 and g 
The form factor for the resonance, Fji{q'^), is taken as. 



2 

pnn 



2.91. 



(37) 



with Aff = 1.0 GeV. The same type of form factors are also applied to the nucleon pole with 
Atv = 0.8 GeV. The propagators of the exchanged meson, nucleon pole and resonance can 
be written as ItI - 



n/r) 



m 



n/r] 



—I 



kp — nip 



(38) 
(39) 
(40) 



TABLE I: Relevant parameters used in our calculation. The masses, widths and branching ratios 
(BR) are taken from central values of PDG [3] except the BR for N*{1A40) Avr. 



Resonance 


Pole Position BW Width Decay Mode Decay Ratio g'^/'in 


A* (1232)^33 


(1210, 100) 


118 


Ntt 


1.0 


19.54 


iV*(1440)Pii 


(1365, 190) 


300 


Ntt 


0.65 


0.51 








Na 


0.075 


3.20 








Avr 


0.135 


4.30 


A^*(1520)L»i3 


(1510, 110) 


115 


Ntt 


0.6 


1.73 








Np 


0.09 


1.32 








Att 


0.2 


0.01 


A^*(1535)5ii 


(1510, 170) 


150 


iVvr 


0.45 


0.037 








Nri 


0.525 


0.34 








Np 


0.02 


0.15 


A*(1600)P33 


(1600, 300) 


350 


Ntt 


0.175 


1.09 








Avr 


0.55 


59.9 








Ar*(1440)7r 


0.225 


289.1 


A* (1620)^31 


(1600, 118) 


145 


Ntt 


0.25 


0.06 








Np 


0.14 


0.37 








Avr 


0.45 


83.7 


A^*(1650)5ii 


(1655, 165) 


165 


Ntt 


0.775 


0.06 








Nt] 


0.065 


0.026 








Np 


0.08 


0.011 








Att 


0.04 


0.063 


N*{W75)Di5 


(1660, 135) 


150 


Ntt 


0.4 


2.16 








Att 


0.55 


3077.5 


7V*(1680)Fi5 


(1675, 120) 


130 


Ntt 


0.675 


5.53 








Na 


0.125 


4.45 








Np 


0.09 


0.32 








Att 


0.1 


9.39 
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TABLE II: Table U continued. 



Resonance 


Pole Position BW Width Decay Mode Decay Ratio g^/^n 


Ar*(l700)L>i3 


(1680, 100) 100 


iVTT 


0.1 


0.075 






Np 


0.07 


0.043 






Att 


0.04 


0.003 


A* (1700)^33 


(1650, 200) 300 


Ntt 


0.15 


1.02 






Np 


0.125 


0.69 






Att 


0.45 


0.072 


A^*(1710)Pii 


(1720, 230) 100 


Ntt 


0.15 


0.012 






Nr] 


0.062 


0.042 






Na 


0.25 


0.085 






Np 


0.15 


36.1 






Att 


0.275 


0.12 


Af*(1720)Pi3 


(1675, 195) 200 


N-K 


0.15 


0.12 






Nrj 


0.04 


0.28 






Np 


0.775 


190.7 



G„(,) = =M±^, (41) 



GfiD = (43) 



q2 _ 




. -^(^ 




q2 _ 


Ml + iMrTr ■ 






q2 _ 


Ml + iMrTr 



Here Vr is the total width of the corresponding resonance, and G^uiq) and G^ua/^iq) is 
defined as, 

11 2 

G^,u{q) = -gt^u + -l^lu + - 7<.9m) + ^j^Qf^lu, (45) 

G^uap{q) = - ^{gf^agui3 + gf^^g^a) + ^gt^ugafi (46) 

- ^ {l^llagup + lulpgiJia + l^il|3gucx + lulag^ip) , (47) 
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Ml' 



7/. = -l^i + 



(48) 



Because constant width is used in the Breit-Wigner (BW) formula, we adopt the pole posi- 
tions of various resonances for parameters appearing in the propagators. 

The coupling constants appearing in relevant resonances are determined by the empirical 
partial decay width of the resonances taken from PDG [l^, and then we adjust the values 
of cut-off in form factors to fit the data. The relations between the branching ratios of 
the adopted resonances and the corresponding coupling constants squared can be calculated 
straightforwardly with above Lagrangians, and most of them can be found in the appendix 
of Ref. Q. The detailed calculations of QpNR and g^NR from the R -> N p{a) Ntttt 
decay are given in Ref. 23|. The values of coupling constants used in our computation 
are compiled in the Table [U together with the properties of the resonances and the central 
value of branch ratios. It should be noted that we adopt a nearly half of the decay width of 
A^*(1440) Avr in PDG as the recent data favored |6|,y,|24|. 

Then the invariant amplitudes can be obtained straightforwardly by applying the Feyn- 
man rules to Fig. 1. As to the different isospin channels, isospin coefficients are considered. 
We do not include the interference terms among different diagrams because their relative 
phases are not known, and the Valencia model seems to show that such terms are very small. 



III. NUMERICAL RESULTS AND DISCUSSION 

As a starting point, in Fig. [2] we demonstrate our calculated total cross sections of six 
isospin channels compared with the existing data [gI. IqIIiI. 1^ . Our numerical results give 
an overall good reproduction to all six channels. The pre-emission diagrams (see (2), (4), 
(6) in Fig. [1]) tend to be negligibly small, consistent with the Valencia model, so we do not 
include them in our concrete computation. In Fig. [2] we do not show the following negligible 
contributions: double-A^*, A^* ^ Np, N* Ntt, A^*(1520) An, A^*(1650) Att, 
A^*(1675) Att, A^*(1680) An, A^*(1680) Na, A^*(1700) An, A^*(1710) ^ Avr, 
A^*(1710) ^ Na, double- A* (1600), A*(1600) ^ A^vr, double- A* (1620), A*(1620) ^ A^p, 
A*(1620) Nn, double-A*(1700), A*(1700) Np, A*(1700) Nn, and A*(1700) An. 
These terms are minor either because of their small branching ratios of double pion channel 
such as A^*(1535), A^*(1650) and A^*(1700), or belonging to higher partial waves such as 
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A*(1620) and A^*(1675), or lying beyond the considered energies such as A^*(1680), A*(1700), 
A^*(1710) and A^*(1720). It should be mentioned that p-meson exchange is much smaller 
than TT-meson exchange in the available diagrams except for nucleon poles but we still include 
the p-meson exchange in the calculation for the completeness of our model. 

Our results underestimate the data in the close-to-threshold region where the final state 
interactions (FSI) should be relevant. We do not consider the initial state interaction (ISI) 
either, because at present we do not have an unambiguous method at hand to simultaneously 
include the FSI and ISI in our model. The ISI usually has a weak energy dependence, so 



20|. 



adjusting cut-off parameters in the form factors may partly account for it effectively 
We would give some qualitative observations of FSI. 

Next we shall first address the pp — > nrm^-K^ channel because it has negligible A^* con- 
tribution to be more clean. Then we shall discuss other channels and explore the different 
situation at each channel. In the following we assume the same definitions of various dif- 



The Mij 



ferential cross sections as graphically illustrated in the experimental articles 
and Mijk are the invariant mass spectra, and the angular distributions are all defined in the 
overall center of mass system. The Qm is the scattering angle of M, and 5^ is the opening 
angle between i and j particles. The 6*-' (or t?*-' corresponding to Gf defined in Ref. js, 8|) 
is the scattering angle of i in the rest frame of i and j with respect to the beam axis (or the 
sum of momenta of i and j). The values of vertical axis are all arbitrarily normalized. 

A. The channel of pp nmr~^TT~^ 

In this channel, we find that the A Ntt Ntctt term is dominant below 1000 MeV 
while the A — > Att and double-nucleon-pole terms are also important. The A — > Nn — ?► Nmr 
term is not included in the Valencia model Our model seems to overestimate the COSY- 
TOF upper-limit by a factor of around four. The A Avr terms in two models are consistent 
with each other because we use the same coupling constant of ttAA from quark model but 
our double-A term contributes smaller as we use a smaller cut-off parameter in ttNA form 
factor. Between 1000 MeV and 1700 MeV, the contribution of the double-A term is the 
most important one, and the nvr^ invariant mass distribution at 1100 MeV do show a clear 
A peak as can be seen in Fig. [31 We also find that the A — > Ntt — > Ntttt and A — > Avr 
terms are crucial to get the right shape of differential cross sections at 1100 MeV. Though 
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the data is of poor statistics, the vr+vr"'" invariant mass spectrum does not show obvious 
low-mass peak and this is reahzed by the inclusion of the A — )■ Njt Ntttt and A Avr 
in our model. The ^tttt has also a significant improvement compared to the double-A alone. 
These distributions should be very useful to constrain the poorly known coupling constant of 
vrAA. The particular enhancement compared to our model without FSI in the nn invariant 
mass spectrum is probably an indication of strong ^Sonn FSI. 

The contribution from A* (1600) — )■ A^*(1440)7r term has a steep rise and begins to take 
over as the largest one for Tp above 1700 MeV. Besides, at large energies contributions from 
the A* (1600) Att and A* (1620) Avr become significant. So in these energy region of 
pp nniT^TT^, it is a good place to explore the properties of these A* resonances. We would 
like to point out that these behaviors together with the dominance of A — > A^vr Ntttt and 
A Att close to threshold alleviate the isospin problem of the pp — > utitt^tt^ and pp 
ppTT^TT^ channels mentioned at the beginning of our article, because the isospin coefficients 
of these terms in pp — > nriTT^TT^ are bigger than that in pp ppn^TT^ channel and this is 
contrary to the case of double-A. As a result, we get an improvement on the description of 
all isospin channels. 

It should be addressed that it is very useful to pin down the cut-off values in form factors 
of the relevant A and A* contributions using the data of pp — )■ nriTT^TT^ at first, and then 
it makes much easier for us to determined the A^* contributions in other channels. The new 
value of total cross section of pp — > nriTT^TT^ measured at CELSIUS js] is in line with previous 
data and this gives our some confidence on the extracted parameters. Further accurate 
measurements of the pp uutt'^tt^ channel should be very helpful for the improvement of 
the model. 



B. The channel of pp ppvr+vr 

Below 1000 MeV the A^*(1440) Na term is the largest while the A^*(1440) Avr term 
is the second. Of these two terms the cr-meson exchange gives much bigger contribution 
than the vr- meson exchange as depicted in Fig. H] and this shows the importance of isoscalar 
excitation of A^*(1440). The Double-A term is negligible at this low energies as well as the 
A Att term. Contributions from the nucleon pole and A^ Att terms are visible below 
700 MeV. The proton and pion angular distributions in the center of mass system at 650 
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and 680 MeV are trivially isotropic and the model does agree with the measured data j6|. 
So we do not show them here. The differential cross sections at 750, 775, 800 and 895 MeV 
are given in Fig. Fig. El Fig. [7] and Fig. [HI respectively. Our model calculations reproduce 
the published data well and are also compatible to the very preliminary data of CELSIUS 
at 895 MeV ^ as shown in Fig. [HI The role of A^*(1440) Att is evident in invariant mass 
spectrums and the data is fitted better than that including A^*(1440) -t- Na alone. Most 
obviously, the anisotropic shape of -(9^+ is well fitted after including the A^*(1440) — ^ Att 
while A^*(1440) — ?■ Na term is symmetric, so -(9^+ together with and is used to 
determined the ratio of partial decay widths of A^*(1440) Att and A^*(1440) Na js]. 
The results are strongly energy dependent and give a smaller decay width of A^*(1440) — ?► Att 
than that listed in PDG. This is believed to support the breathing mode of A^*(1440) 25 |. 
The FSI is evident in pp invariant mass spectrum, but seems to be much weaker than in the 
case of pp — > nmi^Ti^ channel. 

In the Valencia model the double-A is dominant above 1300 MeV. However, because we 
use smaller cut-off parameter for the vrA^A form factor in order to fit both nmi^Ti^ and 
ppn^TT^ channels, our model shows that A^*(1440) Att begins to take over above 1100 
MeV, and double-A and A^*(1440) Na are also important and comparable. We give the 
differential cross sections at llOOMeV and 1360 MeV in Fig. [9] and Fig. [TOl which can be 
tested by the measured data of CELSIUS 8|. The prominent features are the double hump 
structure in M^+j,- and the upward bend in 5,r+7r- which arise from the A^*(1440) — > Att. 
The Valencia model give very similar results because M^+^- and 5^+7r- are sensitive to the 
appearance of A^*(1440) Att. These seem to somewhat incompatible to the preliminary 
data [7, 8 1 which show the phase space behavior in these two spectrums. The same phenom- 
ena also happen in the channel of pp — )> ppir^n^ at high energies, and we will discuss them 
altogether later. 



C. The channel of pp — > ppn tt 

The iV*(1440) -4- Na term dominates below 1000 MeV and the nucleon pole term also 
gives significant contribution below 800 MeV. The A^*(1440) Att and double-A contribu- 
tions are comparable in this energy region. So it should be cautious to use this channel to 
extract the ratio of the partial decay widths for the decay of A^*(1440). Indeed, the extracted 
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ratios from pp — )■ ppTr^vr" are about one third of those from pp — )■ ppTT^TT~ at the same nom- 
inal mass under the assumption that A^*(1440) dominates in this energy range js, 8|. The 
significant double-A and nucleon pole contributions might account for this discrepancy and 
should be reasonably incorporated in the fit. 

Above 1100 MeV, the double-A term dominates; the A^*(1440) Avr and A^*(1440) ^■ 
Na are also important and give similar contributions. Other contributions are much smaller. 
The most striking feature in this energy region is that a level-off behavior happens in the 
total cross section between 1000 and 1200 MeV, while other channels rise smoothly when in- 
creasing the incident energy. Our model fails to describe this behavior and also overestimates 
the high energy data. It is possible that this shape is caused by the interference of different 
diagrams which are not included in our model, but this would require a peculiar energy de- 
pendence of A^* as shown by the isospin decomposition [g], Q] . Another possible explanation 
is that there maybe exist a steep rise of some kind of contribution when other contributions 
are saturated in this energy region. This happens in the channel of pp — )■ nmx^Ti^ where a 
weak level-off at 1600 MeV is caused by the steep rise of A*(1600) A^*(1440)7r. However, 
this is not the case for the pp pp-K^ir^ channel where A* (1600) — > A^*(1440)7r gives much 
smaller contribution due to the isospin factor. So this problem is left for further clarification. 

In Fig. [HI Fig. [12] and Fig. [131 "we show the differential cross sections of pp — )■ ppir^TT^ at 
beam energies of 775, 895 and 1000 MeV, respectively. The data at 775 MeV are well repro- 
duced and A^*(1440) — > Na is overwhelmingly dominant. Some of the angular distributions 
are sensitive to the presence of the A^*(1440) Avr contribution, and hence can be used to 
determine the partial decay ratios of A^*(1440), although this is somewhat complicated by 
the double-A and nucleon pole contributions as we pointed out earlier. The contribution of 
A^*(1440) — 7- Avr and double-A terms become much clearer at 895 and 1000 MeV, though 
slight discrepancy in invariant mass spectrums at 895 MeV exists between our model and 
the measured data, which remind us that it needs a further improvement in the crossover 
region. The phase space shapes of M^o^o and St,Ot^o begin to appear ever since 895 MeV and 
up to high energies of this channel as shown in Fig. [TH Fig. [15] and Fig. [16] which are the 
differential cross sections of pp ppTc^ir^ at 1100, 1200 and 1300 MeV. However, because 
the influence of the A^*(1440) — > Avr does not decrease much or disappear at these energies, 
our model gives a double hump structure in M^o^o and forward peak in S-j^o^q which are 
contradictory to the CELSIUS data. In order to explain the data, it is required that the 
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A^*(1440) — ?■ Avr shows up at low energies but immediately saturated at about 1000 MeV. 
That is a peculiar energy dependence behavior which does not supported by our model. 
Except for M^^r and ^tttd other spectrums at high energies are well fitted by our model both 
in pp — )■ ppTT^n^ and pp ppTT^-K^. So we would rather conclude that something happens 
in the vrvr system which needs a more thorough investigation as the next step. The vr-vr 
rescattering is found to be negligible at these energies |26l |. 

The effect of FSI is not obvious at low energies compared to our calculated curve but 
enhancement seems to happen at high energies. It is possible that this is related to the 
behavior of tt^tt^ system. 

D. The channels of pp — )■ pmr~^ir^, pn — )■ ppir^ir^ and pn — )• pnir^-K^ 

The A^*(1440) — ^ Na does not present in the pp pmr^iT^ reaction, so the double-A 
term is the most important one in a wide energy range. The A — Avr and A Ntt Nnir 
terms have significant contribution below 800 MeV and also have some contribution at higher 
energies together with the A* (1600) and A* (1620) terms. The agreement with the data is 
very good and the FSI may influence the near-threshold region since our model slightly 
underestimates this part. 

The channel of pn pp7T~7r^ is another reaction where the A^*(1440) — t- Na does not 
contribute. Since the charged meson exchange is allowed in this channel, the iV*(1440) Avr 
term is very important and is of the same order as the double-A term in the whole energy 
region. The contributions from the nucleon pole and A Ntt — > NNtt terms are also quite 
significant near the threshold. Our results reproduce the new bubble chamber data measured 



by KEK 10| very well, but underestimate the old data [14] by about a factor of 5. Since 
the double-A contribution has been well determined by the channel of pp — )■ nnir^TT^, we 
think that the main ambiguity comes from the A^*(1440) An term. If future experiments 
confirm the old data, then the isovector mesons like vr- and p-meson should play more 
important role in the excitation of A^*(1440). On the other hand, the new data of KEK 
support the isoscalar excitation of A^*(1440) which is favored by our model. 

The pn — > pn-K^-K^ channel is interesting because it can shed light on the low mass 
enhancement in M^ttt, known as the ABC effect of double pion production in nuclear fusion 
reactions 0. Below 900 MeV, the A^*(1440) Na is found to be dominant while the 
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double-A and A^*(1440) — t- Avr terms also give some contribution. The nucleon pole and 
— )■ Att terms are also important close to threshold. Above 1000 MeV, the double- 
A term is the most important one and the A^*(1440) gives sizable contribution at high 
energies. The total contribution gives a reasonable description to the new KEK data while 
the underestimation of the data close to threshold may be due to the omission of the pn 
FSI. Similar to the pn — > ppn^Tc^ channel, our model does not favor the old bubble chamber 
data which need large isovector excitation of A^*(1440). Very recently the ABC effect is 
experimentally established in pn — > dir^TT^ at beam energies of 1.03 and 1.35 GeV, and has 
been interpreted as an s-channel double-A resonance [271]. According to the observation of 
our model, the A^*(1440) emerges at these energies so it is necessary to take a further look 
at the mechanism of ABC effect in pn — )■ dir^iT^ reaction. As a matter of fact, it has been 
demonstrated that at beam momentum of 1.46 GeV (corresponding to beam energies of 800 
MeV) where the A^*(1440) is expected to be dominant, the deuteron momentum spectra 
3e reasonably explained by the interference of the A^*(1440) — )■ Na and A^*(1440) — )■ 



can 
Att 
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E. Final State Interaction 

As discussed above, the effect of FSI is anticipated to influence the results close to thresh- 
old where the s-wave is expected to be dominant. Usually the Jost function is used to account 
for the FSI enhancement factor, 

^(^■'" - f^l- («) 

where k is the relative momentum of A^A^ subsystem in the final state. The corresponding 
scattering length and effective range are: 

« + 2 

« = — ^ = -—1^^ 50) 

ap a + p 

with a = -7.82fm and r = 2.79fm for ^5*0 pp interaction, a = 5.42fm and r = 1.76fm for 
^Si isoscalar pn interaction, and a = -18.45fm and r = 2.83fm for ^Sq nn interaction. At 
higher energies, the high partial waves become important and above approximate treatment 
would deteriorate rapidly. Fortunately, the effect of FSI should significantly decrease. So 
we may just ignore it above 1.4 GeV. To investigate the influence of the FSI to the energy 
dependence of cross sections, we assume the Jost function for the FSI and normalize this 
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factor to the unity at the beam energy of 1.4GeV. For the final states with the pn pair, we 
assume it is mainly in the '^5'i isoscalar state. 

Though above prescription is quite rough, the agreement with the data are considerably 
improved. In Fig.[T71 we demonstrate the total cross section below 1.4GeV. The effect of FSI 
can be seen in some of the differential cross sections especially A^A^ spectrums. In Fig. [TS] 
we take pp — )• nmr^Ti^ channel as a typical example. The nn FSI gives a sharp peak in 
the nn spectrums which agree with the pp — )■ nmr^n^ data. However, the data in other 
channels do not favor this sharp peak and this reflects the drawbacks of our formalism. The 
nn FSI slightly improves the fit of the rni^n^ and nmc'^ spectrums but increase the slope of 
SnTT+- FSI has very small influence on other angular distributions. The situations for other 
channels are similar. 

IV. SUMMARY 

In this article, we present a simultaneous analysis of varous isospin channels of double 
pion production in nucleon-nucleon collisions up to 2.2 GeV within an effective Lagrangian 
approach. We study the contributions of various resonances with mass up to 1.72 GeV and 
demonstrate that A^*(1440), A, A* (1600), A* (1620) and nucleon pole constitute the main 
ingredients to reasonably explain the measured data while the contribution of other reso- 
nances are negligible. We suggest that it is necessary to consider other influence such as 
the double-A and nucleon pole contributions when one studies the properties of A^*(1440) 
in the channels of pp — ppir^ii^ and pp ppir^TT^. Our model well describe the measured 
differential cross sections except some vrvr spectra which are left as an open problem. Com- 
pared with the Valencia model, the main differences are: (1) Among 3 major ingredients, 
double-A, A^*(1440) — )■ Avr and A^*(1440) — )• Na terms, considered in the Valencia model, 
our model increases signiflcantly the relative contribution from the A^*(1440) — ?■ Na term 
by reducing the relative branching ratio of A^*(1440) — > Avr and assuming a smaller cut-off 
parameter for the vriVA coupling; (2) In addition, our model introduces signiflcant contribu- 
tions from A — )■ Nrr — Ntttt at energies near threshold and from A* (1600) and A* (1620) at 
energies above 1.5 GeV. Though the model should be improved to reasonably incorporated 
the ISI and FSI, the conclusions reached from our model should be helpful to the future 
experiments to be performed at COSY and HIRFL-CSR as well as further theoretical study 
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on related problems. Our results also give hints to the ABC effect in the pn — > dn^n^ and 
pd ^^He7r°7r° reactions which need to be further explored. 
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FIG. 2: (color online) Total cross sections of NN — > NNtttt. The black solid, red short-dash- 
dotted, blue dashed, orange dotted, green dotted, cyan short-dashed, green dash-dotted, royal 
short-dotted, magenta dash-dot-dotted, pink dotted and bold solid curves correspond to contribu- 
tion from double-A, iV*(1440) Na, N*{U40) ^ Avr, A ^ Avr, A Ntt, A* (1600) Avr, 
A* (1600) A^*(1440)7r, A* (1620) An, nu(30on pole, Avr and the fuh contributions, re- 

spectively. The solid circles and triangles represent the data from Ref.O, 9-111. The open circles 




FIG. 3: Differential cross sections of pp — nmr~^Tr^ at beam energies 1100 MeV. The dashed, 
dotted and solid curves correspond to the phase space, double- A and full model distributions. 
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FIG. 4: The iV*(1440) Att, iV*(1440) Na terms of pp ppTr+Tr~. The dashed, dotted and 
sohd curves correspond to 7r-meson exchange, cr-meson exchange and total contribution. 
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FIG. 5: Differential cross sections of pp ppn'^TT^ at beam energies 750 MeV. Tlie daslied, dotted 
and solid curves correspond to the phase space, A^*(1440) Na and full model distributions, 
respectively. The data are from Ref. 
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FIG. 6: Differential cross sections of pp ppvr+vr at beam energies 775 MeV. The meaning of 
curves is the same as Fig. [5j The data are from Ref. gj. 
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FIG. 7: Differential cross sections of pp ppvr+vr at beam energies 800 MeV. The meaning of 
curves is the same as Fig. [5j The data are from Ref. Qj. 
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FIG. 8: Differential cross sections of pp — )■ pptt~^tt at beam energies 895 MeV. The meaning of 
curves is the same as Fig. [5j The preHminary data are from Ref. 8|. 
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FIG. 9: Differential cross sections ofpp — ?> ppTr~^TT~ at beam energies 1100 MeV. The dashed, dotted 
and sohd curves correspond to the phase space, A^*(1440) —?■ Air and full model distributions, 
respectively. The preliminary data are from Ref. 7|. 
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FIG. 11: Differential cross sections of pp — )■ ppir^Tr^ at beam energies 775 MeV. The meaning of 
curves is the same as Fig. [5j The data are from Ref. [sl. 
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FIG. 12: Differential cross sections of pp — )■ ppir^Tr^ at beam energies 895 MeV. The meaning of 
curves is the same as Fig. [5j The data are from Ref. [sl. 
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FIG. 13: Differential cross sections of pp — >■ ppTr^n^ at beam energies 1000 MeV. The dashed, 
dotted and sohd curves correspond to the phase space, double- A and full model distributions. 
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FIG. 14: Differential cross sections of pp — )> ppn^-K^ at beam energies 1100 MeV. The meaning of 
curves is the same as Fig. [131 The data are from Ref. 
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FIG. 15: Differential cross sections of pp — )> ppn^-K^ at beam energies 1200 MeV. The meaning of 
curves is the same as Fig. [131 The data are from Ref. 
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FIG. 16: Differential cross sections of pp — )> ppn^-K^ at beam energies 1300 MeV. The meaning of 
curves is the same as Fig. [131 The data are from Ref. 
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FIG. 17: (color online) Total cross sections of NN NNtttt. The black solid and red dashed 
curves correspond to the full contributions without and with final state interactions, respectively. 
The data are the same as the Fig. [2] 
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FIG. 18: Differential cross sections of pp — )• nnvr'^vr^ at beam energies 1100 MeV. The dashed, 
dotted and soUd curves correspond to the distributions of the phase space, the full model with and 
without FSI, respectively. The data are from Ref. 0]- 



36 



